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Results on high transverse momentum charged particle emission with respect to the reaction plane 
are presented for Au+Au collisions at ^Sj^„= 200 GeV. Two- and four-particle correlations results 
are presented as well as a comparison of azimuthal correlations in Au-|-Au collisions to those in p-l-p 
at the same energy. Elliptic anisotropy, V2, is found to reach its maximum at pt ~ 3 GeV/c, then 
decrease slowly and remain significant up to pt ~ 7 - 10 GeV/c. Stronger suppression is found in 
the back-to-back high-pt particle correlations for particles emitted out-of-plane compared to those 
emitted in-plane. The centrality dependence of V2 at intermediate pt is compared to simple models 
based on jet quenching. 
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111 high energy heavy-ion colhsions, a high density sys- 
tem consisting of deconfined quarks and gluons is ex- 
pected to be created 0. Energetic partons, resulting 
from initial hard scatterings, are predicted to lose energy 
by induced gluon radiation when propagating through 
the medium This energy loss is expected to depend 
strongly on the color charge density of the created sys- 
tem and the traversed path length of the propagating 
parton. At Brookhaven's Relativistic Heavy Ion Collider 
(RHIC) three different observations related to parton en- 
ergy loss have emerged: strong suppression of the inclu- 
sive hadron production 0,0, [J, strong suppression of the 
back-to-back high-p* jet-like correlation ^^'^ large 

values of the elliptic flow at high pt 0- In non-central 
heavy ion collisions, the geometrical overlap region has 
an almond shape in the transverse plane, with its short 
axis in the reaction plane. Depending on the emission 
azimuthal angle, partons traversing this system, on aver- 
age, experience different path lengths and therefore dif- 
ferent energy loss. It leads to (a) azimuthal anisotropy 
in high pt particle production with respect to the reac- 
tion plane 0, 0| (the second harmonic in the particle 
azimuthal distribution, elliptic flow, is characterized 
by V2 = (cos2((/)— VPi?)).) and (b) to the dependence 
of the high pt 2-particle back-to-back correlations on the 
orientation of the pair. 

In this Letter, using higher-order cumulant analy- 
sis ^3 ^3 9.nd comparing azimuthal correlations mea- 
sured in p -\- p collisions to those in Au-|-Au, we confirm 
strong elliptic flow in mid-central Au-|-Au collisions at 
least up to Pt w 7 GeV/c as qualitatively expected in the 
jet quenching scenario. We further investigate the influ- 
ence of the jet quenching mechanism on high-pt particle 
production with respect to the reaction plane by study- 
ing V2 centrality dependence in the intermediate pt re- 
gion and two-particle azimuthal correlations at different 
angles with respect to the reaction plane. 

The data set consists of about 2 million minimum bias 
and 1.2 million central trigger Au-|-Au events and 11 mil- 
lion p -I- p events at ^s„„= 200 GeV. The measurements 
were made using the Time Projection Chamber 0| of 
the STAR detector which covers pscudorapidity (ry) 
from -1.3 to 1.3. The event centrality in this paper is 
defined by the multiplicity measured at mid-rapidity by 
STAR Tracks used to reconstruct the flow vector, 
or generating function |l3l | in the case of the cumulant 
method, were subject to the sa me q uality cuts as used 
in the y/s^= 130 GeV analysis [lall3; except for the 
low transverse momentum cutoff, which for this analysis 
is 0.15 GeV/c instead of 0.10 GeV/c. 

One of the largest uncertainties in elliptic flow mea- 
surements in nuclear collisions is due to so-called non-flow 
effects - the contribution to the azimuthal correlations 
not related to the reaction plane orientation, such as res- 
onance decays and inter- and intra-jet correlations. The 
importance of these effects can be investigated by com- 
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FIG. 1: (color online) Azimuthal correlations in Au-|-Au col- 
lisions (squares) as a function of centrality (peripheral to cen- 
tral from left to right) compared to minimum bias azimuthal 
correlations in p-l-p collisions (circles). Errors are statistical 
only. 

paring the azimuthal correlations measured in Au+Au 
to those in p + p collisions, where all correlations are 
considered to be of non-flow origin. For such a compari- 
son we evaluate the accumulative correlation of a particle 
from a given pt bin with all other particles in the region 
0.15 < Pt < 2.0 GeV/c and I77I < 1.0 by calculating the 
event average sum: 

(^cos2((?!)pj - (t)i)) ^ Mv2{pt)v2 + {non-flow} (1) 

i 

where is the azimuthal angle of the particle from a 
given Pt bin. The first term in the r.h.s. of Eq.^ rep- 
resents the elliptic flow contribution, where V2{pt) is the 
elliptic flow of particles with a given pt , and V2 is the aver- 
age flow of particles used in the sum; M is the multiplicity 
of particles contributing to the sum. The multiplicity in 
the sum changes with the centrality of the collision, but 
as long as the relative number of particles (per trigger 
particle) involved in non-flow effects does not change, the 
contribution due to these effects is a constant. Compar- 
ing p + p and Au-|- Au collisions one indeed might expect 
some changes in particle correlations: there could be an 
increase in correlations due to a possible increase of jet 
multiplicities in Au-|-Au collisions, or conversely, some 
decrease due to the suppression of high pt back-to-back 
correlations 0. It is difficult to make an accurate esti- 
mate of possible modifications of non-flow effects. The 
fact that at very high pt the p+p results arc very close to 
central Au+Au (shown later by Figure QJ, suggests that 
the modifications are relatively small. 

Figure n shows the azimuthal correlation, Eq. ^ as a 
function of transverse momentum for three different cen- 
trality ranges in Au+Au collisions, as compared to min- 
imum bias p + p collisions. We observe that for the most 
peripheral Au+Au collisions, the azimuthal correlations 
are very similar to minimum bias p + p. In mid-central 



4 



0.25- 



0.2- 



0.15- 



- 4k- 

- ★ 



0.05 



3^ Centrality 20%-60% 



2 4 6 8 10 12 

p, (GeV/c) 

FIG. 2: (color online) V2 of charged particles as a func- 
tion of transverse momentum from the two-particle cumu- 
lant method (triangles) and four-particle cumulant method 
(stars). Open circles show the 2-particle correlation results 
after subtracting the correlations measured in p + p collisions. 
Only statistical errors are shown. 



Au+Au events, the azimuthal correlations are very dif- 
ferent from those in p + p collisions in both magnitude 
and pt-dependencc. Note that at pt = 7 GeV/c, the az- 
imuthal correlations in Au+Au collisions are still many 
standard deviations away from those observed in p + p 
collisions, indicating significant elliptic flow up to these 
momenta. For the most central Au+Au collisions, at 
low-pt the magnitude of the correlations is also different 
from p + p. However, for particles with pt > 5 GcV/c, 
the correlation in Au+Au collisions starts to follow that 
in p + p collisions, suggesting that azimuthal correlations 
become dominated by non-flow effects and that the lat- 
ter arc rather similar in p + p and Au+Au collisions at 
those momenta. The observed non-monotonic central- 
ity dependence of the azimuthal correlation at low and 
moderate pt is strong evidence of elliptic flow, ft is qual- 
itatively different from that expected from intra-jet cor- 
relations among jet fragments |l8l |. 

We also perform a multi-particle cumulant analysis, 
which is much less sensitive to non-flow effects than the 
traditional approach based on two-particle correlations. 
Figure El shows V2 as a function of transverse momen- 
tum for 20%-60% of the total cross-section. The V2 ob- 
tained using the four-particle cumulant method, U2{4}, 
is up to about 20% lower than the value of V2 obtained 
from the two-particle cumulant method. This difference 
could be partially explained by non-flow effects, which are 
greatly suppressed in W2{4}, and by the fluctuation of V2 
itself Flow fluctuations contribute to V2{'2} and 

V2{4:} with different signs. The true V2 lies between t'2{4} 
and approximately the average of i'2{2} and W2{4}. The 
systematic uncertainty is given by these two bounds. For 
the centrality range plotted in Fig. El we find significant 
V2 at least up to pt « 7 GeV/c, well within the region 
where particle production is expected to be dominated 



by parton fragmentation. Two-particle cumulant results 
extend to 12 GeV/c, although at high pt these might 
be dominated by non-flow contributions. Also shown in 
Fig. El by open circles are the 2-particle correlation re- 
sults after subtracting the correlations measured in p + p 
collisions. The comparison of these results to W2{4} in 
the region pt < 4 GeV/c indicates that either the rela- 
tive contribution of non-flow effects is larger in Au+Au 
collisions compared to p + p, or there is a significant fiow 
fiuctuation contribution that would increase the apparent 
V2{2} values and decrease V2{4}. In general we observe 
that V2{pt) reaches a maximum at about 3 GeV/c, con- 
firming results obtained by PHENIX and then slowly 
decreases. 

The energy loss mechanism that leads to azimuthal 
anisotropy at high pt also leads to a distinct feature in 
two particle azimuthal correlations. At high transverse 
momenta, two-particle distributions in the relative az- 
imuthal angle measured in p + p, d+Au, and Au+Au 
collisions at RHIC 0, 13 exhibit a jet-like correlation 
characterized by the peaks at Acj) = (near-side corre- 
lations) and at A(/> = tt (back-to-back). The back-to- 
back peak is found to be strongly suppressed in central 
Au+Au collisions 0. In non-central collisions, the sup- 
pression should depend on the relative orientation of the 
back-to-back pair with respect to the reaction plane. In 
the analysis of the two-particle azimuthal correlations, we 
select trigger particles with 4 < p'"^'® < 6 GeV/c emit- 
ted in the direction of the event plane angle ^2 (in-plane, 
|0trig_^2| ^ 7r/4and |(/)*"8;-*2| > 37r/4) and perpendic- 
ular to it (out-of-plane, 7r/4 < - ^'2! < 37r/4). The 
trigger particles are paired with associated particles sat- 
isfying 2 GeV/c< Pt < The tracks are restricted to 
|?7| < I. To reduce the effect of particles produced within 
a jet on the reaction plane reconstruction, all particles in 
a pseudorapidity region |A?7| < 0.5 around the highest pt 
particle in the event are excluded from the event plane 
determination. In the upper panel of Fig. Owe show the 
azimuthal distributions of associated particles for trig- 
ger particles that are in-plane (squares) and out-of-plane 
(triangles) in midcentral Au+Au collisions. The distribu- 
tions are corrected for the reconstruction efficiency. The 
measured distributions exhibit a strong elliptic flow pat- 
tern similar to that found in the recent analysis at the 

sps in. 

In the presence of elliptic flow the in-plane and out-of- 
plane two-particle azimuthal distributions are given by 
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where v^'^^°'^ and Uj"^ are the elliptic flow of the associ- 
ated and trigger particles, respectively, and (cos 2 A?*) is 
the reaction plane resolution [ll|. For the given central- 
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FIG. 3: (color online) Upper panel: Azimuthal distributions 
of associated particles for trigger particles in-plane (squares) 
and out-of-plane (triangles) for Au+Au collisions at centrality 
20-60%. Open symbols are reflections of solid symbols around 
A(j) = and A(f) = vr. Elliptic flow contribution is shown by 
dashed lines. Lower panel: Distributions after substracting 
elliptic flow, and the corresponding measurement in p -|- p 
collisions (histogram) . 



measured via the reaction plane method. For the esti- 
mate of the systematic uncertainty in the determination 
of the flow contribution, we have varied V2^^°'^ and Wg"^ 
between 0.167 and 0.213 (w2{4} and V2{2} measured in 
the range 2< pt <6 GeV/c). To reduce the systematics, 
a z vertex cut of ±25 cm is applied to p + p events to 
match that in Au + Au events. 

The distributions were fit to Eq.[21in the region 0.75 < 
I A^l < 2.24 rad, with B as the only free parameter, to de- 
termine the amount of background. For the in-plane dis- 
tribution, B = 0.649 ±0.004(stoi.)±0.005(s?/s.), and for 
the out-of-plane, B = 0.638 ± 0.004(siat.) ± 0.002(sys.). 
The systematic errors were estimated from using differ- 
ent ranges of A<j> in the fit. We observe a strong ex- 
cess of two-particle correlations over the correlation pat- 
tern generated by elliptic flow in the region \A(j)\ < 0.75 
for both in-plane and out-of-plane distributions, char- 
acteristic of near-side intra-jet correlations. In the re- 
gion around A(f> = n, we observe an excess for the in- 
plane distribution, but no excess is found for the out- 
of-plane distribution. This is better illustrated in the 
lower panel of Fig. |21 where we show the flow-subtracted 
in-plane and out-of-plane distributions compared to that 
measured in p + p collisions. The level of combinatorial 
background measured in p + p collisions, 0.014 ± 0.001, 
has been subtracted. The near-side jet-like correlations 
measured in Au+Au are similar to those measured in 
p + p collisions. The back-to-back (around A0 = tt) 
correlations measured in Au+Au collisions for in-plane 



trigger particles are suppressed compared to p + p, and 
even more suppressed for the out-of-plane trigger parti- 
cles. For the near angle correlations in the relative az- 
imuthal region |A(/)| < 0.75 rad, the integrals of the az- 
imuthal distributions are 0.078 ± 0.014(siat)t^;^^?(s?/s.) 
in-plane and Omi±O.OU{stat.)tomiisys.) out-of-plane. 
For the back-to-back correlations in the relative az- 
imuthal region |A0 — 7r| < 0.75 rad, the integrals 
arc 0M8±0.0U{stat.)to°Qll{sys.) in-plane and 0.014 ± 
0.0U{stat.)to°°2t{sys.) out-of-plane. Note that the large 
systematic errors in Fig. O (lower panel) , resulting from 
the uncertainty in the subtraction of elliptic flow con- 
tribution, are highly anti-correlated: assuming weaker 
(stronger) elliptic flow results in the upper (lower) sys- 
tematic error bar for dn™/dAc/) and lower (upper) sys- 
tematic error bar for dn°"'/c?A0 distributions. 

A different approach to remove the elliptic flow contri- 
bution to the two-particle distributions is to subtract the 
raw away-side correlations from the near-side correlations 
measured in the same \A(j)\ range (in this case, the elliptic 
flow contribution cancels out). The difference in the cor- 
relation strength, an integral over Acfi region, on the near- 
side {\A(j)\ < 0.75 rad) and the away-side (| A0— 7r| < 0.75 
rad) is measured to be 0.030 ± 0.011{stat.) for the in- 
plane triggers and 0.067±0.011(sfat.) for the out-of-plane 
triggers where the systematic uncertainty due to elliptic 
flow is canceled out, and the remaining systematic uncer- 
tainties are smaller than the statistical errors. Assuming 
similar strength of the near-side correlations in-plane and 
out-of-plane, the observed difference can be attributed to 
the suppression of away-side correlations which depends 
on the reaction plane orientation. 

Although results presented above strongly support the 
jet-quenching scenario qualitatively, the amount of ellip- 
tic flow observed at high pt for collisions at ^s„„= 130 
GcV seems to exceed the values expected in the case 
of complete quenching |2^ . Extreme quenching leads to 
emission of high-pt particles predominantly from the sur- 
face, and in this case V2 would be fully determined by the 
geometry of the collision. This hypothesis can be tested 
by studying the centrality dependence of V2 for high-pt 
particles. 

Figure0|shows V2 in the pt-range of 3-6 GeV/c (where 
1)2 is approximately maximal and constant) versus impact 
parameter. The values of the impact parameters were ob- 
tained using a Monte Carlo Glauber calculation . The 
measured values of V2{4:} are compared to various simple 
models of jet quenching. The upper curve corresponds to 
a complete quenching, in which particles are emitted from 
a hard shell |U |2^ ; this gives the maximum values of V2 
that are possible in a surface emission scenario. A more 
realistic calculation corresponds to a parameterization of 
jet energy loss in a static medium where the absorption 
coefflcient is set to match the suppression of the inclu- 
sive hadron yields 0. The density distributions of the 
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FIG. 4: (color online) ?;2 at 3 < pt < 6 GeV/c versus impact 
parameter, b, compared to models of particle emission by a 
static source (see text). 

static medium are modeled using a step function (follow- 
ing [i^) and a more realistic Woods-Saxon distribution 
(following [2^). The corresponding V2 values are shown 
as the upper and lower band, respectively. The lower and 
upper boundaries of bands correspond to an absorption 
that gives a suppression factor of 3 and 5 0, respec- 
tively, in central collisions. Over the whole centrality 
range, the measured V2 values are much larger compared 
to calculations. Taking into account that this measure- 
ment is dominated by the lower pt side (3 GeV/c), the 
quark coalescence mechanism [28| might be responsible 
for the difference, but no quantitative explanation for the 
observed large elliptic flow exists at the moment. 

In summary, we have shown that the charged parti- 
cle elliptic anisotropy in midcentral Au-|-Au collisions at 
^s„„ =200 GeV extends to large transverse momenta, at 
least up to ~ 7 GeV/c, as expected in a jet quenching 
scenario. By performing multi-particle correlation analy- 
sis and comparing the azimuthal correlations in Au-t-Au 
collisions to those in p+p, we find the contribution of the 
effects not associated with the reaction plane orientation 
is relatively small in midcentral events but could be sig- 
nificant in peripheral and central collisions. We report 
stronger suppression of the back-to-back high pt correla- 
tions for out-of-plane triggers compared to in-plane trig- 
gers, again consistent with a jet quenching picture. V2 
integrated from moderate to high pt, approximately in 
the region where it reaches a maximum, clearly exceeds 
the limits set for elliptic flow due to a simple jet quench- 
ing mechanism, and still waits for quantitative theoretical 
explanation. 

We thank the RHIC Operations Group and RCF at 
BNL, and the NERSC Center at LBNL for their sup- 
port. This work was supported in part by the HENP 
Divisions of the Office of Science of the U.S. DOE; the 
U.S. NSF; the BMBF of Germany; IN2P3, RA, RPL, 
and EMN of France; EPSRC of the United Kingdom; 
FAPESP of Brazil; the Russian Ministry of Science and 
Technology; the Ministry of Education and the NNSFC 



[1 
[2: 

[s: 
[4: 
[s: 
[6: 
[7: 

[8 

[9: 

[10: 
[11 

[12 

[13 

[14 

[15 

[16 

[17 

[18 

[19 
[20 

[21 

[22; 

[23' 
[24; 



[25' 
[26; 

[27 
[28; 



See, for example, J.W. Harris and B. Miiller, 
Annu. Rev. Nucl. Part. Sci. 46, 71 (1996). 
M. Gyulassy and M. Plumer, Phys. Rev. Lett. B 243, 432 
(1990); X.N. Wang and M. Gyulassy, Phys. Rev. D 44, 
3501 (1991); Phys. Rev. Lett. 68 1480 (1992); R. Baier, 

D. Schiff and B. G. Zakharov, Ann. Rev. Nucl. Part. Sci. 
50, 37 (2000). 

K. Adcox et al. [PHENIX Collaboration], Phys. Rev. 
Lett. 88, 022301 (2002). 

C. Adler et al. [STAR Collaboration], Phys. Rev. Lett. 

89, 202301 (2002). 

J. Adams et al. [STAR Collaboration], Phys. Rev. Lett. 
91, 172302 (2003). 

C. Adler et al. [STAR Collaboration], Phys. Rev. Lett. 

90, 082302 (2003). 

J. Adams et al. [STAR Collaboration], Phys. Rev. Lett. 

91, 072304 (2003). 

C. Adler et al. [STAR Collaboration], Phys. Rev. Lett. 
90, 032301 (2003). 

R. J. M. SneUings, A. M. Poskanzer and S. A. Voloshin, 
nucl-Gx/9904003 

X. N. Wang, Phys. Rev. C 63, 054902 (2001). 

A. M. Poskanzer and S. A. Voloshin, Phys. Rev. C 58, 
1671 (1998). 

N. Borghini, P. M. Dinh, and J.-Y. OUitrauh, Phys. Rev. 
C 63, 054906 (2001). 

N. Borghini, P. M. Dinh, and J.-Y. OUitrauh, Phys. Rev. 
C 64, 054901 (2001). 

M. Anderson et al, Nucl. Instrum. Meth. A 499, 659 
(2003). 

K. H. Ackermann et al. [STAR Collaboration], Nucl. In- 
strum. Meth. A 499, 624 (2003). 

K. H. Ackermann et al. [STAR Collaboration], Phys. 
Rev. Lett. 86, 402 (2001). 

C. Adler et al. [STAR Collaboration], Phys. Rev. C 66, 
034904 (2002). 

Y. V. Kovchegov and K. L. Tuchin, Nucl. Phys. A 708 

413c (2002); Nucl. Phys. A 717 249c (2003). 

M. Miller and R. SneUings, nucl-ex/0312008 

S.S. Adler et al. [PHENIX Collaboration], Phys. Rev. 

Lett. 91, 182301 (2003). 

G. Agakichiev et al. [CERES/NA45 Collaboration], 
Phys. Rev. Lett. 92, 032301 (2004). 

J. Bielcikova, S. Esumi, K. Filimonov, S. Voloshin, and 
J. P. Wurm, Phys. Rev. C 69, 021901(R) (2004). 

E. V. Shuryak, Phys. Rev. C 66, 027902 (2002). 

B. B. Back et al, Phys. Rev. C 65, 031901(R) (2002); 
K. Adcox et al, Phys. Rev. Lett 86, 3500 (2001); I. G. 
Bearden et al, Phys. Lett. B 523, 227 (2001). 

S. A. Voloshin, Nucl. Phys. A 715, 379c (2003). 
Xin-Nian Wang, private communication. Calculation is 
based on the framework of Phys. Lett. B 595, 165 (2004). 
A. Drees, H. Feng, and J. Jia, nucl-th/0310044 

D. Molnar and S. Voloshin, Pliys. Rev. Lett. 91, 092301 
(2003). 



